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1.  Summary 


The  ground  vehicle  suspension  system  has  a  very  strong  influence  on  the  handling,  stability,  ride, 
and  durability  characteristics  of  a  vehicle  as  it  negotiates  different  types  of  terrains.  Therefore, 
research  advancements  in  suspension  and  control  system  technologies  are  needed  to  innovate 
better  systems,  which  are  compact,  self-powered,  easily  adaptable,  and  tunable  to  different 
vehicle  weight  configurations  and  different  terrain  types.  Military  ground  vehicles  with  high  a 
center  of  gravity  and  large  variations  in  gross  weights  can  greatly  benefit  by  using  advanced 
suspension  technologies  with  active,  adaptive,  preview-sensor-based  control  algorithms.  The 
ultimate  objective  of  this  research  work  is  to  develop  a  preview-sensor-based  nonlinear  control 
algorithm  that  can  effectively  work  with  semi-active  dampers,  such  as  the  magnetorheological 
(MR)  fluid  damper.  This  research  study  attempts  to  develop  advanced  suspension  technology 
with  active,  adaptive  control  algorithms  using  dynamics  simulation  tools  and  then  validate  them 
experimentally  to  demonstrate  the  desired  effectiveness  by  reducing  vehicle  body  acceleration  by 
10%  or  increasing  off-road  speed  by  10%.  This  report  discusses  and  summarizes  the  research 
activities  conducted  by  the  U.S.  Army  Research  Laboratory,  Vehicle  Technology  Directorate 
(ARL-VTD).  The  main  tasks  of  this  research,  to  be  accomplished  over  a  span  of  two  years,  are 
as  follows: 

•  Develop  and  adapt  generalized  predictive  control  algorithms  and  extend  them  to  include 
nonlinear  system  identification  for  use  with  MR  dampers. 

•  Characterize  performance  of  developed  control  algorithms  and  select  the  optimum 
algorithm  based  on  simulation  using  a  quarter-car  suspension  system  and  multibody 
vehicle  dynamics  models. 

•  Determine  effectiveness  of  control  algorithms  based  on  the  ability  to  control  vibration 
caused  by  terrain  disturbances,  as  well  as  overall  vehicle  handling  and  ride  performance. 

•  Leverage  ARL  broad  agency  agreement  (BAA)  research  with  the  University  of  Maryland 
to  develop  an  advanced  prototype  MR  damper  that  can  verify  the  effectiveness  of  control 
algorithms  using  a  quarter-car  representative  2-degree-of-freedom  (2-DOF)  Advanced 
Suspension  Test  Rig  (ASTR)  being  currently  installed  at  ARL-VTD ’s  laboratory. 

•  Build  test  setup  and  characterize  baseline  and  advanced  MR  dampers  on  quarter-car 
representative  ASTR.  Conduct  verification  and  validation  of  control  algorithms  to  assess 
their  vibration  control  effectiveness. 

This  report  specifically  covers  the  tasks  conducted  during  the  first  year  of  research  to  develop  an 
advanced  semi-active  suspension  system  and  adaptive  control  algorithms  for  improving  vehicle 
vibration  control  and  ride  quality  of  military  ground  vehicles. 


1 


2.  Introduction 


Passive  suspension  systems  for  Army  wheeled  vehicles  are  optimized  for  passenger  comfort  and 
road/terrain  handling  qualities.  These  two  goals  often  lead  to  conflicting  design  requirements. 
With  a  view  toward  designing  suspension  systems  that  would  improve  both  performance  metrics, 
ARL-VTD  has  underway  a  research  program  aimed  at  developing  active  suspension  systems  for 
implementation  on  the  next  generation  of  Army  vehicles.  The  overall  objective  is  to  reduce 
vibration-induced  fatigue  as  well  as  improve  target  aiming  precision  in  theater. 

There  are  basically  two  main  groups  of  active  suspension  systems:  semi-  and  fully  active.  Semi¬ 
active  suspension  systems  vary  spring  and  damper  properties  to  add  or  remove  energy  from  the 
system  according  to  some  applied  control  algorithm.  Usually,  this  is  achieved  through  variable 
damping  components,  such  as  variable  orifice  valves  or  MR  fluids.  Fully  active  suspension 
systems  use  force  actuators  to  add  and  remove  energy  from  the  system.  This  function  can  be 
achieved  using  hydraulic,  pneumatic,  or  electromechanical  actuators.  Current  research  at  VTD  is 
addressing  both  types  of  systems. 

The  generalized  predictive  control  (GPC)  algorithm  is  a  linear,  time-invariant,  multi-input/multi- 
output  method  that  uses  an  AutoRegressive  exogenous  (ARX)  model  to  represent  the  input- 
output  behavior  of  the  system  and  to  design  the  controller  (1-4).  Three  variants  of  the 
techniques  have  been  evaluated.  In  the  first  variant,  no  explicit  account  is  taken  of  the 
disturbances  acting  on  the  system  during  system  identification  and  controller  design  (2,  3).  This 
version  is  called  GPC  With  Implicit  Disturbance  in  this  report.  The  second  and  third  variants 
take  explicit  account  of  measurable  disturbances.  The  first  of  these  two  (called  GPC  With 
Explicit  Disturbance  herein)  uses  only  the  past  values  of  the  measured  disturbances,  while  the 
latter  (called  GPC  With  Preview  in  this  report)  additionally  uses  measurements  (or  predictions) 
of  future  disturbances  (1,4). 

This  report  discusses  the  control  methodology  of  GPC  algorithms,  vehicle  dynamics  models  and 
co-simulation  with  control  algorithms,  extension  of  GPC  to  enable  nonlinear  capability, 
advanced  MR  damper  development,  and  2-DOF  suspension  test  rig  assembly  and  setup.  The 
report  finishes  with  some  concluding  remarks  on  the  results  of  the  completed  research  activities 
and  some  comments  on  the  nature  of  some  planned  future  work. 
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3.  Vehicle  Dynamics  Model 


TruckSim  is  a  widely  used,  commercially  available  software  tool  used  to  simulate  and  analyze 
the  dynamic  behavior  of  wheeled  vehicles.  It  is  employed  in  conjunction  with  Matlab/Simulink 
to  simulate  the  effects  of  control  algorithms  on  a  full-vehicle  model.  The  TruckSim  model  takes 
the  full  vehicle  and  terrain  into  consideration,  including  driver  controls  and  steering,  road 
geometry,  suspensions,  tires,  and  power  train.  Component  parameters  are  input  into  TruckSim 
via  constants,  linear  coefficients,  nonlinear  tables,  and  algebraic  formulae. 

Full  three-dimensional  road  geometry  was  specified  at  0.1 -m  increments  in  this  research  work. 
The  steering  system  captures  the  steer  of  each  wheel  due  to  the  steering  system  geometry  and 
compliances.  The  suspension  system  models  capture  the  full  nonlinear  kinematical  behaviors 
and  compliances  of  the  front  and  rear  suspensions.  These  models  include  data  for  springs, 
dampers,  and  jounce/rebound  stops.  The  tire  model  is  an  internal  table-based,  single-point 
contact  model.  Nonlinear  tables  using  actual  measured  tire  data  represent  vertical  force,  lateral 
force,  longitudinal  force,  aligning  moment,  and  overturning  moment  as  functions  of  deflection, 
slip,  load,  and  camber.  The  power-train  model  includes  engine  torque,  torque  convertor 
characteristics,  and  transmission  gear  ratios  and  efficiencies.  A  representative  TruckSim  tactical 
vehicle  model  (modified  High-Mobility  Multipurpose  Wheeled  Vehicle  [HMMWV]  model)  that 
was  developed  in  collaboration  with  the  U.S.  Army  Materiel  Systems  Analysis  Activity  (5)  is 
shown  in  figure  1.  Figure  2  contains  a  graphical  layout  of  the  suspension  as  modeled  in 
SuspensionSim  that  works  with  TruckSim  software. 


Figure  1 .  TruckSim-modified  HMMWV  model  on 
Perryman3. 
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Figure  2.  Representative  HMMWV  suspension 
model. 

For  this  research,  we  modified  some  of  the  basic  HMMWV  vehicle  parameters  to  study  the 
performance  of  GPC  algorithms  on  various  terrain  profiles.  The  vehicle  mass  was  doubled.  The 
suspension  springs  were  stiffened,  and  the  front  damping  coefficient  was  reduced.  This  is  shown 
in  table  1. 


Table  1.  Vehicle  parameters  for  modified  HMMWV  model. 


Parameter 

Basic 

HMMWV 

Modified 

HMMWV 

Sprung  mass 

2210  kg 

5000  kg 

Front  spring 

150  N/mm 

250  N/mm 

Front  damper 

30  kNs/m 

15  kNs/m 

Rear  spring 

150  N/mm 

400  N/mm 

Rear  damper 

30  kNs/m 

30  kNs/m 

4.  Control  Algorithm  Methodology 


The  steps  in  developing  a  GPC-based  controller  are  system  identification  and  derivation  of  the 
control  law.  System  identification  is  based  on  the  derivation  of  a  multistep  output  prediction 
equation  using  the  observer  Markov  parameters  (OMPs)  that  comprise  the  ARX  model  used  to 
characterize  the  system.  The  control  law  then  follows  from  the  minimization  of  a  quadratic 
performance  equation  involving  the  predicted  response  and  the  control  command.  A  summary  of 
the  key  steps  of  the  process,  adapted  from  the  more  extensive  development  present  in  Kvaternik 
et  al.  (4),  is  given  next.  Since  GPC  is  a  multi-input/multi-output  controller,  the  control 
command, {ut},  measured  response, {yt},  and  external  disturbance,  {dt},  can  be  vectors  whose 
elements  correspond  to  the  values  of  a  specific  actuator  or  sensor  at  time  t. 
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4.1  GPC  With  Implicit  Disturbance 


The  relationship  between  the  input  and  output  time  histories  of  a  linear  time-invariant  multi- 
input/multi-output  system  can  be  described  by  an  ARX  model  of  the  form 

y(t)  =  ajy(t  -  1)  +  a2y(t  -  2)  +  ...  +  avy(t  -  p ) 

+  / 30u(t )  +  /?i u(t  —  1)  +  ...  +  Ppu(t  -  p)  ^ 

This  equation  says  that  the  output,  y(t),  at  the  current  time  step,  t,  may  be  estimated  by  using  p 
sets  of  the  previous  input  and  output  measurements,  u(t-l),  u(t-p),  y(t-l),  ...,y(t-p),  and  the 

current  input,  u(t).  The  integer  p  denotes  the  order  of  the  ARX  model.  The  coefficient  matrices 
a;  and  Pi  are  called  observer  Markov  parameters  and  are  the  quantities  computed  by  the 
identification  process  (<5). 

To  start  the  identification  process,  the  system  is  excited  with  band-limited  white  noise.  These 
independent  random  excitations  are  applied  concurrently  to  all  of  the  rc  control  inputs.  The  rc 
control  inputs  and  the  m  response  outputs  are  recorded  for  L  samples.  The  resulting  input  and 
output  time  histories,  u  and  y,  are  then  used  to  form  the  data  matrices  y  and  V. 


y  =  YV, 


where 


y=  [Fo  -  Vl- i] 

mxL 


u0  u± 

U2 

Up 

UL- 1 

u0 

Mi 

•  Up- 1 

UL- 2 

yo 

Fi 

-  yP- 1 

yi- 2 

u„ 

-  Up— 2 

UL- 3 

yo 

•  yv- 2 

UL- 3 

u0 

•  uL_p_1 

. 

yo 

y  L—p—i- 

[rc  +  (rc  +  m)p]  x  L 


(2) 

(3) 


(4) 


Throughout  this  report,  the  dimensions  of  the  matrices  and  vectors  are  noted  in  blue  below  key 
equations.  The  order  of  the  ARX  model,  p,  and  the  number  of  samples,  L,  are  specified  by  the 
user.  Some  guidelines  for  their  selection  are  given  in  Juang  (6).  The  matrix  Y,  containing  the 
OMPs,  a;  and  Pi,  follows  from  equation  2  by  inversion  and  has  the  following  form: 

F  =  [Pa  Pi  p2  «2  ■■■  Pp  «p]  , 

mx[rc  +  (rc  +  m)p] 
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The  one-step-ahead  output  prediction  equation  given  by  equation  1  is  the  starting  point  for 
deriving  the  multistep  output  prediction  equation  that  is  needed  for  deriving  the  GPC  controller. 
This  equation  is  obtained  by  replacing  t  by  t+j  in  equation  1,  and  letting  /  range  over  the  set  of 
values  j=l,  2,  ...  hp-l,  where  hp  is  the  prediction  horizon  (the  number  of  time  steps  for  which  the 
future  responses  are  predicted).  The  resulting  equations  can  be  assembled  into  a  multistep  output 
prediction  equation  having  the  following  form: 

jy/J  =  x{uhJ  +  «(yp}  +  p{up} .  (6) 

The  coefficient  matrices  x,  p,  and  a  are  formed  from  combinations  of  the  OMP.  The  quantity 
yh  is  the  vector  containing  the  future  predicted  responses,  whereas  uh  is  the  vector  containing 
the  (unknown)  future  control  commands.  The  integer  hc  is  the  control  horizon  (the  number  of 
time  steps  over  which  future  control  is  assumed  to  act).  The  quantities  up  and  yp  are  vectors 
containing  the  previous  p  sets  of  control  commands,  and  measured  responses,  respectively. 


yt+ o 


[yt+hp- ij 


mhp  x  1 


Po 

Po 


Po 


q-1 


PS 

Po 


L/?oP  1 


PS 

PS 


q-2 


q-1 


Po 


hp  2 


Po 


Po 


hp-hc 


mhp  x  rrh, 


c“-c 


mhp  x  mp 


Pi 

Pp 

Pi 

Pi 

-  Pi 

pr 

pr1 

-  pr1 

pi 

pq 

-  pqv 

.pTx 

pr1 

■ 

mhp  x  rcp 


rchc  x  1 


ai 

a2 

Up  ~ 

a\ 

a 

q-1 

Up 

ai 

a\ 

cc q 

Lip 

hp- 1 
<*1 

hp—1 

a2 

hp  —  1 

CC 

up  J 

mp  x  1 


(Ut-l] 


[u 


t-p) 


rcp  x  1 


(7) 
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(8) 

(9) 


where  the  coefficients  are  calculated  using  the  following  relations: 

Pi  =  pr  +  ■ 

q  q-1  .  q— 1 

c(p_1  =  a*  +  ap_t . 

Pi-i  =  pv  +  «rvP-i  •  do) 

In  this  formulation  the  size  of  the  ARX  model,  p,  is  set  equal  to  the  number  of  future  predicted 
responses,  hp,  (prediction  horizon)  and  future  commands,  hc,  (control  horizon).  Generally,  these 
three  parameters  can  differ  from  one  another  (i,  2). 

The  optimal  control  law  is  obtained  by  minimizing  an  objective  function.  To  accomplish  this, 
one  defines  a  quadratic  objective  function  as  follows: 

J  =  yThpQyhp  +  uThcRuhc ,  (li) 

where  Q  is  a  symmetric  positive  semi-definite  matrix  that  assigns  weights  to  the  predicted 
responses,  R  is  a  symmetric  positive  definite  matrix  that  assigns  weights  on  the  future  control 
commands,  and  the  superscript  T  indicates  the  vector  or  matrix  transpose. 

Substituting  the  expanded  equation  for  y  yields 

/  =  Mw/d  +  abP}  +  p{up\]t  QWiuhc}  +  «{yP]  +  P{wP}]  +  uThRuhc .  (i2) 

The  minimum  is  obtained  by  taking  the  partial  derivative  of  the  function  with  respect  to  future 
commands, [uhc],  and  setting  it  equal  to  zero. 

0  =  [tKJ  +  «frp}  +  P{«p}]r<?[xKc}  +  “bp)  +  P{“P}]  +  J^cukRuK  ■  (13) 

(  Ut+  0  \ 

0  =  2tt Q[x{uhc]  +  a{yp}  +  p{up}]  +  2r\  \ .  (14) 

(Wt+ftJ 

Solving  for  control  commands  yields  the  following  optimal  future  commands: 

[uh J  =  -(xrQx  +  R)-1xr(?[a{yp}  +  p(up}] .  (15) 

Defining  the  following  control  law  matrices, 

Ac  =  —(ttQt  +  R)_1xrQa  (16) 

and 

Bc  =  -(xr(?x  +  R)-1xrQ/?,  (17) 

yields  control  laws  that  calculate  the  next  hc  sets  of  commands: 

{«hj  =  Acbp}  +  Bc{up}  ,  (18) 

where  uhc  corresponds  to  the  control  commands  hc  time  steps  into  the  future. 
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However,  only  the  first  set  of  control  commands  will  be  applied  to  the  system,  thus  only  the  first 
rc  rows  of  {uhc }  are  needed.  The  remaining  rows  are  discarded. 


(Vt- i) 

(ut- 1) 

{«t+0} 

=  «c 

;  +Pc  ; 

yt-p) 

\Ut-p) 

rcxl 

mp  x  1 

rcp  x  1 

rc  x  mp  rc  x  rcp 


(19) 


4.2  GPC  With  Explicit  Disturbance 

Generalized  Predictive  Controller  With  Explicit  Disturbance,  as  discussed  in  Juang  and  Eure  (7) 
and  Kvaternik  et  al.  (4),  is  similar  to  GPC  With  Implicit  Disturbance,  except  that  now  the 
external  disturbances  are  explicitly  measured  and  included  in  the  model.  Time  histories  of 
control  commands,  measurement  responses,  and  now  disturbance  measurements  are  used  to 
generate  the  ARX  model  for  the  system.  The  V  matrix  with  disturbance  information  is  shown 
next. 


d  o 

d'L 

dp 

dL-i 

u0 

tip 

-  U-L-l 

do 

dp-i 

d[j—2 

u0  • 

tip  _  i 

UL- 2 

To  •' 

••  yv- 1 

-  Tl-2 

u0 

V-L-p-l 

_ 

y0 

VL-p-l- 

[rd  +rc  +  (rd  +rc  +  m)p]  x  L 

The  new  matrix  F,  containing  the  OMP,  has  the  following  form: 

Y  =  [So  /?o  Sj  /?x  a-y  52  Pz  az  ■■■  Sp  /?p  ap] 

mx[rd  +  rc  +  (rd  +  rc  +  m)p ] 


(20) 


(21) 


These  OMP  values  are  used  to  predict  future  responses  based  on  past  responses  (including 
external  disturbances)  and  control  commands. 

jy/J  =  +  a(yp}  +  p{up}  +  5{dp} . 


f  yt+o 

[yt+hp 


'  Ut+ 0  ) 

i  PV 

;  1 

■  +  a]  : 

Mt+hc) 

1  (yt-P 

(Ut-U 


+p 


u 


t-p 


+  8 


(22) 

(23) 


As  before,  the  next  control  commands  are  calculated  by  minimizing  the  objective  function  and 
taking  the  partial  derivative  with  respect  to  control  commands;  the  rows  corresponding  to  future 
control  commands  are  dropped.  Using  this  extension,  we  can  use  the  measured  disturbance 
information. 
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Tr  X  1 


(24) 


fVt-i 

(“t+oj  =  ac  I  : 

{yt-P) 


+  Pc(  ;  )  +  6c 


U 


t-pj 


mp  x  1  rcp  x  1 

rc  x  mp  rc  x  rcp 


rdp  x  1 
rc  x  rdp 


4.3  GPC  With  Preview  Control 

Preview  Control  or  Look-Ahead  Control  measures  the  external  disturbance  before  it  affects  the 
vehicle  and  incorporates  this  information  into  the  control  algorithms.  Preview  Control  was  first 
published  in  Bender  (7)  and  is  further  developed  in  Van  der  Aa  (8)  and  El  Madany  et  al.  (9)  for 
linear  full-state  feedback  control. 


Generalized  Predictive  Controller  With  Preview  Control  is  a  direct  extension  of  GPC  With 
Explicit  Disturbance.  If  it  is  possible  to  generate  a  model  of  the  explicit  disturbance  to  predict 
future  disturbances  (4),  then  this  future  information  can  be  used  in  conjunction  with  the  control 
command,  external  disturbance,  and  measured  response  time  histories  to  generate  the  future 
control  commands.  This  research,  however,  assumes  that  the  future  disturbance  can  be  directly 
measured  using  vehicle-based  sensors.  Again,  time  histories  of  controls  commands, 
measurement  responses,  and  disturbance  measurements  are  used  to  generate  the  ARX  model  for 
the  system.  The  matrix  Y  is  identical  to  equation  21. 


These  time  histories  are  used  to  predict  future  responses  based  on  past  responses,  past  and  future 
disturbances,  and  past  and  future  control  commands.  In  this  formulation,  we  assume  that  the 
current  and  future  disturbance  can  be  measured.  The  matrix  y  is  the  component  of  the  ARX 
model  that  corresponds  to  the  current  and  future  disturbances. 


ryt+o 

Wt+hp 


'  ut+0 

.ut+hc 


ut-i 

;  +  p< 

(yt-pj 

Mt-p . 

(25) 


As  before,  the  future  control  commands  are  calculated  by  minimizing  the  objective  function  and 
taking  the  partial  derivative  with  respect  to  control  commands;  only  the  rows  corresponding  to 
the  next  immediate  control  commands  are  retained.  Using  this  extension,  we  can  use  both  the 
future  and  the  previous  measured  disturbance  information. 


{Mf+o}  =  Yc 


■yt-v 

Ut- 

(dt- 1 

+  «c 

jt-p. 

+  Pc 

Mt-p, 

+  SC  i 
Mt-p 

(26) 


The  difference  in  control  command  between  GPC  With  Explicit  Disturbance  and  GPC  With 
Preview  Control  is  solely  due  to  the  presence  of  the  future  disturbance  term.  The  control 
matrices,  ac,  Pc,  and  5c  are  the  same  for  both  control  laws. 
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5.  Road  Profiles 


Three  road  profiles  were  used  to  evaluate  the  algorithms:  Sine,  Curb,  and  Perryman3.  In  each 
case,  the  road  profiles  were  provided  as  an  array  of  ground  distances  and  vertical  heights. 

5.1  Sine  Road  Profile 

This  road  profile  was  defined  by  a  sine  wave  with  a  half-amplitude  of  0.10  m  and  a  wavelength 
of  10  m.  The  intent  was  to  provide  a  single-frequency  disturbance  that  relates  to  vehicle  speed — 
in  this  case,  1  m/s  (3.6  km/h)  of  speed  corresponds  to  0.1  Hz  excitation.  There  is  a  30-m  flat 
stretch  to  allow  the  model  to  reach  steady  state  after  startup.  Figure  3  shows  a  portion  of  this 
type  of  road  profile. 


Road  Distance  -  m 


Figure  3.  Sine  road  profile. 

5.2  Curb  Road  Profile 

This  road  profile  is  defined  by  a  0.1 0-m  step  rise  followed  by  10  m  of  level  surface,  and  then  a 
0.10-m  step  drop,  as  shown  in  figure  4.  The  intent  was  to  provide  a  step  input  disturbance  into 
the  system.  There  is  a  30-m  flat  stretch  to  let  the  initial  oscillations  die  out  and  the  vehicle  reach 
steady  state. 
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Figure  4.  Curb  road  profile. 


5.3  Perryman3  Road  Profile 

Perryman3  is  an  off-road  test  track  located  at  the  Aberdeen  Proving  Ground.  This  road  profile, 
shown  in  figure  5,  is  intended  to  evaluate  the  algorithms  against  a  challenging  “real-world” 
course. 


Figure  5.  Perryman3  road  profile. 


6.  Simulation  Results:  Linear  GPC  Variants 


Simulations  were  run  until  the  vehicle  traveled  150  m.  The  open-loop  model  was  run  in  parallel 
with  closed-loop  GPC  With  Implicit  Disturbance,  GPC  With  Explicit  Disturbance,  and  GPC 
With  Preview  Control  for  performance  comparison  studies.  All  measurements  assumed  no  noise 
or  offset  bias.  Initial  conditions  for  the  suspension  system  state  variables  were  set  to  zero. 
Vehicle  speed  was  held  at  a  constant  value  for  the  entire  run. 
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For  this  work  the  system  identification  for  each  algorithm  was  precomputed.  Band-limited  white 
noise  was  used  to  excite  the  four  actuators.  The  vehicle  was  run  over  the  Perryman3  terrain 
profile  at  40  km/h.  This  provided  broadband  random  disturbance  inputs  out  to  50  Hz  with  large 
peaks  from  0.5  to  3  Hz. 

The  following  objective  functions  were  used  as  metrics  in  several  comparisons  of  algorithm 
performance: 

jResp  =  lowpassiyffy)  ,  (27) 

Actuator  Force  =  low  pass(VtFti  )  ,  (28) 


and 


Actuator  Power  =  low  pass (^[(power^power  )  ,  (29) 

where  “low  pass”  denotes  a  fourth-order  digital  low-pass  Butterworth  filter  with  a  cutoff 
frequency  of  1  Hz. 

Using  the  cosine  double  angle  formula,  we  see  that  the  square  of  a  signal  will  have  a  steady 
offset  component  and  a  high-frequency  component  that  is  twice  the  frequency  of  the  original 
signal.  The  low-pass  filter  removes  the  oscillatory  component,  leaving  only  the  steady  offset, 
which  is  related  to  the  amplitude  of  the  fast  Fourier  transform  of  the  original  signal.  Because  the 
cutoff  frequency  is  higher  than  the  lowest  frequency  of  excitation  in  the  road  profiles,  there  is  a 
small  low-frequency  oscillation  in  the  objective  function,  J.  The  low-frequency  oscillation  is 
removed  by  averaging  the  objective  function  over  the  entire  run  time. 


For  each  simulation,  the  averaged  values  of  J  were  plotted  for  a  series  of  vehicle  speeds  up  to 
80  km/h.  The  closed-loop  responses  were  normalized  by  the  corresponding  values  of  the  open- 
loop  responses  at  each  speed. 

Response  Improvement  =  1  -  .  (30) 

jResp,open  loop 


This  normalization  allows  a  more  meaningful  comparison  of  the  algorithms  despite  the  variation 
in  actuator  force  or  power. 


Force  Normalized  _  Response  Improvementaig0rithm 
Eff  ectiveness  Actuator  Force aig0rithm 

Power  Normalized  Response  Improvement  algorithm 

Eff  ectiveness  Actuator  Power  algorithm 


(31) 

(32) 


Simulation  results  based  on  the  modified  HMMWV  parameters  are  presented  in  figures  6-12. 
These  figures  were  generated  using  the  following  algorithm  parameters,  as  shown  in  table  2. 
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Figure  6.  Effect  of  control  weights  on  GPC  with  preview  control  on  Sine  road  profile. 


Figure  7.  Effect  of  GPC  algorithms  on  Sine  road  profile. 
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Figure  8.  Time  history  of  control  algorithms  for  light  vehicle  parameters  on  Sine  course  at  50  km/h. 


Figure  9.  Effect  of  GPC  algorithms  for  modified  HMMWV  on  Perryman3  road  profile. 
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Figure  10.  Time  history  of  GPC  algorithms  for  modified  HMMWV  parameters  on  Perryman3  course  at  50  km/h. 


Figure  11.  Effect  of  GPC  algorithms  on  Curb  road  profile. 


15 


33  3.6  3.7  3.8  3.9  4.0 

Time  -  s 


33  3.6  3.7  3.8  3.9  4.0 

Time  -  s 


Figure  12.  Time  history  of  GPC  algorithms  for  light  vehicle  parameters  on  Curb  course  at  72  km/h. 


Table  2.  Algorithm  tuning  parameters  used  for  modified  HMMWV  vehicle. 


Parameter 

GPC  With 
Implicit 
Disturbance 

GPC  With 
Explicit 
Disturbance 

GPC  With 
Preview 
Control 

Actuators,  rc 

4 

Disturbances,  rd 

4 

Sensors,  m 

3 

Q 

50  (vertical  acceleration),  10  (pitch),  1  (roll) 

R 

range  {  le-7  to  5e-6  } 

ARX  model  size,  p 

8 

Prediction  horizon,  hp 

8 

Control  horizon,  hc 

8 

Preview  horizon 

— 

— 

8  (0.08s) 

ID  sample  size,  L 

400  points 

Disturbance  for  ID 

Perryman  3  at  40  km/h 

Actuator  ID  signal 

White  noise 

Actuator  limit 

1.5e5  N 
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6.1  Effects  of  Control  Weighting 

The  modified  HMMWV  model  was  run  with  GPC  With  Preview  Control  over  the  Sine  road 
profile  from  0  to  80  km/h. 

The  objective  function,  actuator  force,  and  actuator  power  are  averaged  over  a  150-m  simulation. 
The  speed  range  corresponds  to  a  sinusoidal  excitation  of  0  to  2.2  Hz.  This  is  repeated  several 
times  for  a  range  of  control  weights  (R  matrix  values);  the  response  weights  (Q  matrix  values) 
were  held  constant.  This  effect  of  the  control  weighting  matrix  is  shown  in  figure  6. 

An  order  of  magnitude  decrease  in  the  control  weight  matrix  significantly  improves  the  response 
and  increases  the  required  force  and  power.  The  force-normalized  effectiveness  shows  that  the 
curves  for  each  weight  coalesce  in  the  same  regions;  the  exception  of  R  =  0.1  e-6  occurs  because 
the  forces  required  are  at  the  limit  for  the  actuator.  This  suggests  that  as  long  as  the  commanded 
forces  are  within  the  capability  of  the  actuator,  the  level  of  the  response  is  linear  with  control 
loads.  However,  when  comparing  the  actuator  power,  the  power-normalized  effectiveness 
decreases  as  the  control  weights  decrease.  This  shows  that  the  actuator  power  costs  increase 
much  faster  than  the  responses  decrease.  Actuator  power,  then,  becomes  the  driving  factor  for 
choosing  control  weights. 

6.2  Comparison  of  Linear  GPC  Algorithms 

A  short  comparison  of  the  three  control  laws  operating  on  the  three  road  courses  (disturbance 
inputs)  is  presented  in  the  following  discussion.  The  control  weight  for  GPC  With  Explicit 
Disturbance  was  R  =  1.0e-7  over  the  entire  speed  range.  To  simplify  analysis,  the  control 
weights  for  the  other  two  algorithms  were  chosen  such  that  the  magnitude  of  the  measured 
responses  matched  GPC  With  Explicit  Disturbance  at  each  speed.  Again,  the  objective  function 
and  the  actuator  force  are  averaged  over  the  150-m  run  for  sinusoidal  disturbances  from  0  to 
80  km/h  (corresponding  to  0  to  2.2  Hz).  Figure  7  compares  the  three  algorithms  vs.  vehicle  speed. 

The  open-loop  response  reaches  a  peak  at  50  km/h,  representing  the  vehicle’s  pitch  mode 
resonance  frequency.  The  closed-loop  algorithms  reduce  the  peak  by  60%,  effectively 
attenuating  the  natural  frequency.  The  normalized  data  shows  that  GPC  With  Explicit 
Disturbance  and  GPC  With  Preview  Control  have  similar  performance,  requiring  the  same 
actuator  force  and  power.  However,  GPC  With  Implicit  requires  significantly  greater  actuator 
force  and  power — almost  twice  the  actuator  power  at  80  km/h.  The  addition  of  road  disturbance 
information  to  the  GPC  algorithm  greatly  reduces  the  actuator  loads  for  the  closed-loop  system. 

Figure  8  is  a  time  history  of  the  measured  responses  and  actuator  forces  for  each  algorithm  as  the 
vehicle  maneuvers  over  the  Sine  road  profile  at  50  km/h,  which  represents  a  single-frequency 
disturbance  of  1 .4  Hz. 
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All  three  versions  of  the  GPC  algorithm  reduce  the  chassis  vertical  acceleration  and  pitch. 
Chassis  roll  increases  because  of  the  very  low  response  weight  assigned  to  chassis  roll.  GPC 
With  Implicit  Disturbance  performs  worse  on  controlling  chassis  vertical  acceleration  but 
slightly  better  at  reducing  chassis  pitch.  Since  the  chassis  pitch  is  the  dominant  term  in  the 
objective  function,  similar  levels  result  for  the  response  objective  function.  It  is  clear  that  GPC 
With  Implicit  Disturbance  requires  much  higher  actuator  loads  than  the  other  two  algorithms. 

The  same  simulations  were  run  using  the  Perryman3  road  profile;  the  algorithm  performance  is 
shown  in  figure  9.  The  GPC  algorithms  reduce  the  overall  response  by  more  than  25%  from  the 
highest  open-loop  response.  Much  like  the  performance  on  the  Sine  road  profile,  GPC  With 
Explicit  Disturbance  and  GPC  With  Preview  Control  have  almost  identical  response  and  actuator 
objective  functions.  Once  again,  GPC  With  Implicit  Disturbance  requires  significantly  higher 
actuator  forces  and  loads  to  match  the  response  objective  function. 

The  GPC  algorithms  reduce  the  peaks  of  the  chassis  pitch  compared  to  open-loop  cases.  This 
dominant  sensor  makes  up  for  the  slight  increases  in  closed-loop  roll  and  vertical  acceleration. 
Without  measured  knowledge  of  the  disturbance,  GPC  With  Implicit  Disturbance  must  wait  until 
the  sign  of  the  responses  changes  before  applying  a  countering  force.  This  change  results  in  the 
large,  but  brief,  forces  as  the  algorithm  catches  up  to  the  changing  slope  of  the  disturbance.  The 
algorithms  with  disturbance  information  have  a  much  quicker  response,  thus  reducing  the  peaks. 
The  time  history  of  the  vehicle  on  Perryman3  is  shown  in  figure  10. 

The  performance  of  the  algorithms  for  the  Curb  road  profile  is  shown  in  figure  1 1 .  Again,  all 
three  GPC  algorithms  reduce  the  response  objective  function  by  as  much  as  14%  at  the  highest 
speed.  However,  this  time  there  is  a  surprising  difference  between  the  actuator  force  and  power 
of  the  algorithms.  GPC  With  Explicit  Disturbance  requires  the  highest  actuator  loads,  double  the 
force  and  power  as  required  for  GPC  With  Implicit  Disturbance.  GPC  With  Preview  Control 
requires  the  least  power  of  the  three  variants — less  than  200  W  above  50  km/h,  as  shown  in 
figure  1 1 . 

Previously,  GPC  With  Explicit  Disturbance  was  nearly  indistinguishable  from  GPC  With 
Preview  Control.  Figure  12  provides  a  closer  look  at  the  actuator  forces  experienced  as  the 
vehicle  encounters  the  step  height  increase  of  the  road.  When  the  GPC  With  Explicit 
Disturbance  encounters  a  step  increase  in  the  disturbance  height  (basically,  a  very  steep  slope), 
this  steep  slope  produces  a  high  control  command.  The  steep  slope  is  in  the  disturbance  history  for 
eight  time  steps. 

On  the  other  hand,  GPC  With  Preview  Control  issues  small  commands  (100  N)  for  the  eight  time 
steps  prior  to  the  start  of  the  curb.  This  reduces  the  peak  force  required  once  the  tire  contacts  the 
curb,  resulting  in  the  lowest  peak  force.  Additionally,  after  the  tire  contacts  the  curb,  the  future 
disturbance  contains  eight  time  steps  of  flat  data  to  further  temper  the  control  commands. 
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As  before,  GPC  With  Implicit  Disturbance  is  unaware  of  the  disturbance  and  can  only  respond 
after  the  measured  responses  are  affected  by  the  curb.  Thus,  it  commands  high  peaks  but  not  in 
response  to  an  infinite  slope. 


7.  Simulation  With  Nonlinear  GPC 


The  research  focus  was  to  demonstrate  nonlinear  system  identification  and  control  of  a  quarter- 
car  suspension  system.  The  requirement  for  nonlinearity  is  driven  by  the  use  of  semi-active  MR 
dampers,  which  are  inherently  nonlinear.  The  algorithm  that  was  considered  and  developed  for 
this  task  is  called  bilinear  generalized  predictive  control  (bGPC).  The  bGPC  algorithm  is  an 
extension  of  the  previous  work  completed  at  ARL-VTD  on  GPC  with  an  implicit  disturbance 
model.  The  bGPC  algorithm  includes  two  additional  terms  in  comparison  to  GPC — namely, 
future  commands  coupled  with  future  responses  and  the  past  commands  coupled  with  past 
responses.  The  NASA  report  by  Jer-Nan  Juang  (10)  describes  the  formulation  of  this  type  of 
model  predictive  control  algorithm.  The  results  of  this  system  ID  and  controller  (bGPC)  were 
obtained  using  simulation  with  a  quarter-car  suspension  model  with  MR  damper.  The  vehicle 
parameters  from  a  generic  vehicle  model  from  previously  conducted  research  (11)  at  ARL-VTD 
were  used  in  a  quarter-car  suspension  model,  as  shown  in  figure  13.  The  predictions  from 
dynamic  simulations  of  the  bGPC  algorithm  for  a  Sine  wave  road  profile  are  shown  in  figure  14. 
The  effects  on  chassis  acceleration  and  percentage  improvement  between  open-loop  and  bGPC 
closed-loop  cases  are  shown  in  figure  14.  Also,  a  full  vehicle  multibody  dynamics  model 
developed  using  TruckSim  was  used  to  analyze  the  effectiveness  of  the  bGPC  algorithm  through 
cosimulation.  Figure  15  shows  the  performance  of  the  bGPC  algorithm  on  a  modified  HMMWV 
model  (TruckSim  model)  for  the  Perryman3  road  profile.  The  simulation  results  (figure  15) 
show  -11%  improvement  for  the  most  challenging  Perryman3  terrain  profile. 

A  different  type  of  nonlinear  control  algorithm,  called  single  experiment-multiple  pulses,  was 
also  considered  and  developed  using  Matlab/Simulink.  This  algorithm  uses  a  series  of  impulses 
to  excite  the  command  inputs  and  generate  a  system  ID  from  those  responses.  For  this  algorithm 
development,  a  simplified  linear  quadratic  regulator  was  used  and  adapted  to  implement  closed- 
loop  control.  However,  this  control  algorithm  was  unstable  in  system  ID  simulation  and  hence 
was  dropped  from  the  research  effort. 
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Figure  14.  Comparison  of  bGPC  algorithm  with  quarter-car  model  for  Sine  road  profile. 
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Figure  15.  Performance  of  bGPC  vs.  linear  GPC  and  open-loop  cases  (modified  HMMWV  model  with 
Perryman3  profile). 


8.  Joint  Research  With  University  of  Maryland  on  Advanced  MR  Damper 


Semi-active  shock  and  vibration  isolation  systems  using  MR  energy  absorbers  (MREAs)  require 
minimization  of  the  field-off  damping  force  at  high  speed  (disabling  additional  stiffness 
generated  by  magneto-rheological  fluid  effects).  This  is  because  the  viscous  damping  force  for 
high  shaft  speed  becomes  excessive  and  will  transmit  the  terrain-induced  vibration  to  the  vehicle 
body.  This  effect  indicates  that  the  controllable  dynamic  force  range,  defined  as  the  ratio  of  the 
field-on  damping  force  to  the  field-off  damping  force,  is  dramatically  reduced.  In  addition,  fail¬ 
safe  MREA  performance,  if  power  were  lost,  is  important  to  shock  and  vibration  isolation 
systems.  A  key  design  goal  is  to  minimize  the  field-off  damping  force  while  maximizing  MREA 
dynamic  force  and  also  maintaining  fail-safe  performance.  This  ARL  BAA  research  has  come 
up  with  the  principle  of  a  bidirectional,  controllable  MREA  that  can  produce  large  damping  force 
and  dynamic  force  range,  as  well  as  excellent  fail-safe  performance.  The  structural  design  of  this 
device  consists  of  electromagnetic  coils  and  permanent  magnets  alternatively  arranged  in  the 
inner  cylinder  of  the  device,  as  shown  in  figure  16.  The  permanent  magnetic  fields  generated  by 
the  permanent  magnets  can  be  changed  or  even  totally  cancelled  by  applying  negative  current  to 
the  electromagnetic  coils.  A  particular  advantage  of  this  magnetic  circuit  is  that  a  fail-safe 
damping  force  in  the  event  of  lost  power  is  provided  by  the  permanent  magnets.  The  damping 
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force  of  the  MREA  can  be  increased  (or  decreased)  when  applying  positive  (or  negative)  current 
to  the  electromagnetic  coils  to  strengthen  (or  weaken)  the  magnetic  field  strength  in  the  magnetic 
flux  path,  which  realizes  bidirectional  control  of  the  MREA  damping  force.  The  University  of 
Maryland  fabricated  a  baseline  prototype  of  this  MREA  device.  This  MR  damper  will  be  used  in 
the  ASTR  at  ARL-VTD  for  carrying  out  verification  and  validation  of  developed  control 
algorithms. 


Figure  16.  Bidirectional  controllable 
inner  bypass  MR  damper 
with  fail-safe  capability. 


9.  Advanced  Suspension  Test  Rig  (ASTR) 


A  2-DOF  ASTR  has  been  fabricated  to  provide  a  basis  for  validating  control  laws  that  are 
developed.  The  test  rig  has  variable  sprung  and  unsprung  masses  of  up  to  250  kg,  representing  a 
1000-kg  ground  vehicle.  The  tire  and  suspension  components  can  be  reconfigured  using 
commercial  off-the-shelf  components  or  custom-made  prototype  components  as  needed.  This 
test  rig  will  be  modified  during  FY  2013  to  adapt  it  for  use  on  the  MTS  load  frame  testing 
machine  that  ARL-VTD  recently  purchased.  This  load  frame  machine  will  be  used  to  input  road 
profile  disturbances  into  the  suspension  system  and  study  the  effectiveness  of  control  algorithms. 
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The  suspension  damper,  shown  in  the  figure  17,  is  replaceable  with  advanced  active  or  semi¬ 
active  dampers  depending  on  our  research  interests.  dSPACE  hardware  equipment  has  been 
purchased  as  well  to  conduct  hardware-in-the-loop  experiments  for  validating  control  laws.  The 
control  algorithm  verification  and  validation  tests  are  planned  to  be  conducted  during  FY  2013. 
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Figure  17.  ASTR  for  quarter-car  representation. 


10.  Conclusions 


This  section  summarizes  the  results  of  numerical  simulations  directed  at  evaluating  the 
effectiveness  of  GPC  for  the  active  control  of  chassis  suspension  systems  in  Army  wheeled 
vehicles.  GPC  is  a  linear  time-invariant,  multi-input/multi-output  predictive  control  method  that 
uses  an  ARX  model  to  describe  the  input-output  relationship  characterizing  the  system.  The 
coefficient  matrices  of  the  ARX  equation  are  determined  using  system  identification  techniques 
and  are  then  used  to  form  the  matrices  comprising  the  GPC  control  law.  Key  equations  in  the 
formulation  of  the  method  were  summarized. 

Three  variants  of  the  basic  GPC  linear  algorithm  were  investigated.  The  first  assumes  that  the 
disturbances  are  unknown  and  computes  the  control  law,  taking  implicit  account  of  the 
disturbances  acting  on  the  system  (GPC  With  Implicit  Disturbance).  The  second  assumes  that 
the  past  values  of  the  disturbances  can  be  measured  and  uses  that  information  in  deriving  the 
control  law  (GPC  With  External  Disturbance).  The  third  variant  assumes  that  the  disturbances 
ahead  of  the  vehicle  can  be  measured  and  additionally  uses  that  data  when  deriving  the  control 
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law  (GPC  With  Preview  Control).  Closed-loop  behavior  of  the  system  was  simulated  using  a 
full-vehicle  nonlinear  model  of  a  two-axle  tactical  vehicle  moving  over  three  types  of  road 
profile  disturbances  for  each  of  the  three  GPC  variants  and  compared  to  the  open-loop  behavior 
of  the  system.  These  comparisons  show  significant  reductions  in  the  chassis  motions  of  the 
model. 

Simulation  has  shown  that  for  a  given  algorithm,  the  control  weights  reduce  sensor  objective 
function  at  the  same  rate  that  they  increase  actuator  objective  function.  Thus,  the  control  weight 
can  be  directly  used  to  scale  the  total  force  available  to  a  controller.  GPC  With  Implicit 
Disturbance  consistently  requires  higher  actuator  forces  and  more  power  than  GPC  With  Preview 
Control  or  GPC  With  Explicit  Disturbance  when  given  a  tonal  or  broadband  disturbance.  GPC 
With  Preview  Control  performs  nearly  identically  to  baseline  GPC  With  Explicit  Disturbance; 
however,  the  preview  control  algorithm  does  drastically  reduce  the  initial  pulse  due  to  step 
inputs.  Despite  the  presence  of  disturbance  information,  GPC  With  Explicit  Disturbance 
performs  poorly  on  step  road  profile,  even  worse  than  having  no  disturbance  information  at  all. 

In  summary,  for  the  three  road  course  profiles  investigated,  GPC  With  Preview  Control  performs 
better  overall  than  the  other  GPC  variants  for  reducing  vehicle-chassis  acceleration  and  pitch. 

The  simulation  results  show  that  linear  GPC  With  Preview  Control  is  capable  of  reducing  the 
vehicle-chassis  acceleration  by  as  much  as  25%  on  average  for  the  most  challenging  terrain 
profile  course  simulated  with  active  suspension  system.  However,  the  nonlinear  version — 
namely,  bGPC — shows  -11%  improvement  for  the  most  challenging  Perryman3  terrain  profile. 
The  differences  in  the  performance  of  the  various  algorithms  studied  in  this  research  effort  are 
currently  being  examined  further  to  understand  the  main  causes  for  the  observed  performance 
differences.  In  the  future,  these  GPC  algorithms  will  be  tested  using  a  2-DOF  suspension  test  rig 
at  the  Army  Vehicle  Research  Laboratory  to  verify  and  validate  these  algorithms. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2-DOF 

AMSAA 

ARL 

ARX 

ASTR 

BAA 

bGPC 

d 

F 

GPC 

HMMWV 

hp,  hc 

Hz 

km/h 

L 

m 

m,  rc,  rd 

MR 

MREA 

N 

OMP 

P 

Q,R 

s 


2  degrees  of  freedom 

U.S.  Army  Materiel  Systems  Analysis  Activity 
U.S.  Army  Research  Laboratory 
AutoRegressive  exogenous  input 
Advanced  Suspension  Test  Rig 
broad  agency  agreement 
bilinear  generalized  predictive  control 
disturbance  inputs 
control  actuator  force 
generalized  predictive  control 
High-Mobility  Multipurpose  Wheeled  Vehicle 
prediction  and  control  horizons 
hertz 

kilometers  per  hour 

number  of  samples  used  for  system  identification 
meter 

number  of  response  outputs,  control  inputs,  disturbance  inputs 

magnetorheological 

magnetorheological  energy  absorber 

newton 

observer  Markov  parameters 
order  of  the  ARX  model 
response  penalty,  control  penalty 
seconds 
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t 

u 

u(t),  y(t),  d(t) 
u(t-j)  ,y(t-j),  d(t-j) 
u(t+j)  ,y(t+j),  d(t+j) 
V 

VTD 

XS,  Xy 

%S>  %U>  Zl  >  %2 
XSlXu 

y 

7 

zl’z2 

on.  Pi  5i 
t^c?  Pc?  Yc?  5C 
x,  a,  p,  y,  5 


index  of  current  time  step 
control  inputs 

vectors  of  current  inputs,  outputs,  and  disturbances  at  time  step  t 

vectors  of  past  inputs,  outputs,  and  disturbances  at  time  step  (t-j) 

vectors  of  future  inputs,  outputs,  and  disturbances  at  time  step  ( t+j ) 

system  identification  data  matrix 

Vehicle  Technology  Directorate 

absolute  coordinates  of  sprung  and  unsprung  masses 

velocity  of  xs,  xv,  z1(  z2 

acceleration  of  xs,  xv 

response  outputs 

matrix  of  observer  Markov  parameters 
compression  of  the  suspension  system  and  the  tire 
observer  Markov  parameters 
control  law  gain  matrices 

coefficient  matrices  in  multi-step  output  prediction  equation 
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